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THE  EFFECT  OF  UNIFORM  FLOW  ON  THE  DYNAMICS  AND  ACOUSTICS 
OF  FORCE-EXCITED  INFINITE  PLATES 


BY 

D  J  ATKINS 


.Summary 

^Formulae  are  obtained  for  the  time-harmonic  steady-state 
plate  displacement,  acoustic  pressure  and  particle  velocities. 

The  £ar-£ield  pressure  is  obtained  from  a  stationary  phase 
approximation,  and  two  separate  formulations  are  used  to  define  a 
time-averaged  acoustic  intensity  vector.  The  stability  of  a  line- 
excited  plate  is  examined  using  both  a  steady-state  approach  and 
a  causal  approach  which  is  able  to  distinguish  between  absolute 
and  convective  instability .  Some  numerical  results  are  presented 
for  a  steel  plate  vibrating  in  water.  The  far-field  pressure  is 
not  significantly  affect <?d  by  flow  speeds  of  a  practical 
magnitude.  The  very  low-frequency  acoustic  intensity  vectors  at 
the  plate  surface  are  affected  considerably,  but  at  higher 
frequencies  there  is  no  significant  effect.  At  low  flow  speeds  a 
convective  instability  is  found  and  at  higher  flow  speeds  there 
is  an  absolute  instability . 


AMTE  (Teddington) 
Queen's  Road 

TEDDINGTON  Middx  TW11  OLN 


62  pages  December  1982 

19  figures 

C 

Copyr ight 

Controller  HMSO  London 
1982 


-1- 


1.  Introduction 


7 


2.  Point-excited  infinite  plate  9 

3.  Line-excited  infinite  plate  11 

4.  Definition  of  the  acoustic  intensity  vector  12 

5.  Far-field  radiated  sound  pressure  14 

(a)  Point-force  excitation  14 

(b)  Line-force  excitation  lb 

6.  Discussion  of  stability  16 

(a)  Steady-state  approach  16 

(b)  Causal  approach  18 

7.  Numerical  results  23 

(a)  Dispersion  plots  23 

(b)  Far-field  pressure  computations  using  steady-state 

approach  25 

(c)  Intensity  vector  plots  using  steady-state  approach  25 

(d)  Numerical  determination  of  stability  characteristics 

using  causal  approach  26 

8.  Discussion  of  results  28 

9.  Conclusions  and  suggestions  for  further  research  29 

Figures  1-19  35 

Tables  I  and  II  53 

Appendix  A.  Derivation  of  the  far-field  sound  pressure 

of  a  point-excited  plate  55 


Appendix  B.  Far-field  sound  pressure  of  a  point-excited 
plate  expressed  in  traditional  spherical 
co-ordinates  59 


-3- 


c 


CP 
Co 
D 
E 
f 

Fo 

h 

i 
-* 

I 

I x ' ly  > ^z 
lm[  ] 
k 

kP 

M 

N 

Nx , Ny , Nz 

P 

P 

R 

Re[  ] 

t 

U 


vx ' vy » vz 
vx,vy,vz 

w(x,y, t) 
W(x,y) 


W(x) 


W(x,t) 


W 

-x* 

W 


speed  of  sound  in  fluid 

phase  speed  of  free  waves  in  plate-fluid  system 
phase  speed  of  free  waves  along  plate  in  vacuo 
bending  stiffness  of  plate  [«Eh3/{12 ( 1-v2 ) } ] 
plate  Young's  modulus 
frequency  in  Hz 

point  force  or  line  force  per  unit  length 

plate  thickness 

V(-l) 

acoustic  intensity  vector 

components  of  acoustic  intensity  vector 

imaginary  part 

wavenumber  in  fluid  (~w/ c) 

free  plate  wave  number  {»(pQhw2/D) 1^A } 

Mach  number  (-*U/c) 

energy  flux  vector 

components  of  energy  flux  vector 

instantaneous  value  of  acoustic  pressure 

acoustic  pressure  amplitude 

radial  distance  from  excitation 

real  part 

time 

flow  speed 

instantaneous  values  of  acoustic  particle  velocity 

acoustic  particle  velocity  amplitudes 

instantaneous  value  of  plate  displacement 

plate  displacement  amplitude  (point-force 
excitation) 

plate  displacement  amplitude  (line-force 
excitation) 

plate  displacement  amplitude  (line-force 
excitation,  causal  approach) 

Pourier  transform  of  W  or  W 

Laplace  transform  of  W 


■*.j  _ 


x,y,  2 
Z(w,a,£) 
Z(u,a) 
a,0 


•V 


y(a,/3) 

0(  ) 

C 

V 

9 


I 


V 

p 

a 

ir 

W,CJ0 

* 


rectangular  Cartesian  co-ordinates 

dispersion  function  for  a  point-excited  plate 

dispersion  function  for  a  line-excited  plate 

Fourier  transform  wave  numbers  (a  only  for 
line-excited  plates) 

Fourier  transform  wavenumber  of  acoustic 
pressure  and  acoustic  velocity  potential 

delta  function 

dimensionless  wavenumber 

plate  loss  factor 

angle  in  cylindrical  and  spherical  co-ordinate 
system 

fluid-loading  parameter  (»p/ps h) 

Poisson's  ratio  for  plate 

fluid  density 

plate  material  density 

Imlw]=o  denotes  line  of  Laplace  integration 
instantaneous  value  of  acoustic  velocity  potential 
amplitude  of  acoustic  velocity  potential 
angle  in  spherical  co-ordinate  system 
radian  frequency  (*2rrf) 

(as  superscript)  denotes  complex  conjugate 


Much  theoretical  work  has  been  carried  out  on  the  dynamics 
of  point-  and  line-excited  infinite  fluid-loaded  plates,  with  the 
fluid  assumed  to  be  stationary  [1-7].  In  practice  however,  either 
the  plate  or  the  fluid  may  be  in  motion  and  hence  an 
understanding  of  the  effect  of  flow  on  the  properties  of  such 
dynamical  systems  is  required.  Here,  expressions  are  derived  for 
the  plate  response  and  radiated  sound  pressure  for  point-  and 
line-excited  infinite  plates  in  the  presence  of  a  uniform 
inviscid  flow.  The  expressions  are  obtained  as  one-  and  two- 
dimensional  Fourier  integrals  over  a  wavenumber  space,  by  solving 
the  steady-state  equations  of  motion  using  Fourier  transforms. 

There  has  been  much  debate  in  the  literature  over  the  past 
twenty  years  on  a  theoretical  definition  of  acoustic  intensity  in 
the  presence  of  flow,  with  particular  application  to  duct 
acoustics.  There  has  also  been  a  recent  interest  in  the 
measurement  of  acoustic  intensity  rather  than  just  the  sound 
pressure  level  as  such  measurements  indicate  the  direction  of 
acoustic  energy  flow.  At  AMTE,  Spicer  [8]  has  numerically 
computed  and  plotted  acoustic  intensity  vectors  in  the  near  field 
of  a  line-excited  plate  without  flow.  It  is  difficult  to  attach  a 
meaning  to  the  term  "acoustic  intensity"  in  situations  where 
there  is  mean  flow.  However,  two  definitions  of  the  acoustic 
intensity  vector  have  been  proposed  and  discussed  in  the 
literature  [9-12],  both  of  which  satisfy  a  conservation  relation, 
i.e.  the  net  outflow  of  acoustic  energy  from  any  volume  equals 
the  rate  of  decrease  of  the  acoustic  energy  contained  within  that 
volume.  In  the  present  study,  the  work  of  Spicer  has  been 
extended  to  include  the  effect  of  uniform  flow.  Intensity  vectors 
are  evaluated  numerically  using  both  definitions  and  the  results 
ar  e  compar  ed . 

In  the  absence  of  flow,  the  only  source  of  energy  in  the 
plate/fluid  system  is  the  exciting  force.  Hence  the  net  acoustic 
energy  flow  is  outwards  from  this  force.  Without  plate  damping, 
free  waves  must  either  propagate  or  decay  away  from  the 
excitation.  A  moving  fluid  possesses  energy  and  hence  energy 
transfer  between  the  plate  and  the  fluid  or  vice  versa  may  give 
rise  to  instabilities.  The  general  problem  of  the  effect  of  a 
flexible  boundary  on  hydrodynamic  stability  has  been  studied  by 
Benjamin  [13,14],  Landahl  [15],  Miles  [16]  and  others.  These 
workers  found  theoretically  that  there  are  three  different  types 
of  instability  possible,  denoted  by  them  as  classes  A,  B  and  C. 
Class  A  disturbances  are  the  Tollmein-Schlicht ing  waves 
encountered  in  stability  treatments  of  flow  over  a  rigid  surface; 
these  disturbances  are  modified  by  the  presence  of  a  flexible 
boundary  and  destabilized  by  damping.  Class  B  disturbances  are 
free  surface  waves  which  can  propagate  along  the  boundary  if  the 
fluid  is  absent  but  which  are  modified  by  the  reaction  of  the 
disturbed  flow.  Class  C  disturbances  are  analogous  to  the  Kelvin- 
Helmholtz  instability  associated  with  vortex  sheets  and  these 


disturbances  occur  when  conservative  hydrodynamic  forces  cause  a 
unidirectional  transfer  of  energy  to  the  boundary.  Gupta  et  al. 

[17]  in  a  laboratory  experiment  observed  class  A  and  class  B 
instabilities  on  an  air-water  interface.  Class  A  waves  could  be 
excited  at  low  air  speeds  by  a  ribbon  vibrating  in  the  air  and 
the  class  B  waves  at  slightly  higher  speeds  by  a  ribbon  vibrating 
in  the  water.  There  was  some  evidence  of  class  C  instability  at 
still  higher  air  speeds. 

An  alternative  approach  to  stability  was  adopted  by  Melcher 

[18]  in  consideration  of  wave  propagation  in  flows  with  electric 
and  magnetic  fields.  He  applied  the  principle  of  causality,  i.e. 
that  the  system  remains  undisturbed  until  at  some  instant  an 
excitation  is  suddenly  applied,  and  the  system  response  is  then 
calculated  after  a  long  time  interval.  There  are  three 
possibilities:  first  that  the  response  at  all  points  of  the 
system  consists  of  components  represented  by  propagating  or 
decaying  waves,  secondly  that  the  response  at  all  points  may 
increase  exponentially  with  time,  and  thirdly  that  the  response 
contains  a  propagating  wave  which  increases  in  amplitude  with 
distance  although  the  response  at  a  fixed  point  does  not  amplify 
with  time.  These  three  possibilities  are  usually  termed  "stable", 
"absolutely  unstable"  and  "convectively  unstable"  respectively. 

Brazier-Smith  &  Scott  [19],  using  a  causal  approach,  have 
examined  theoretically  the  stability  of  uniform  incompressible 
flow  over  a  line-excited  undamped  infinite  plate.  They  deduced  a 
critical  flow  speed  above  which  the  system  is  absolutely  unstable 
and  below  which  the  downstream  response  exhibits  a  convective 
instability.  In  the  present  work,  similar  results  have  been 
obtained  numerically  for  compressible  flow,  and  the  effect  of 
plate  damping  on  stability  is  considered.  The  causal  approach  and 
the  Benjamin-Landahl  approach  to  stability  are  compared  and  shown 
to  give  similar  results  at  low  flow  speeds,  but  different  results 
at  high  flow  speeds. 
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2.  POINT-EXCITED  INFINITE  PLATE 


A  uniform  thin  plate  lies  in  the  plane  z=0  and  it  is 
subjected  to  a  time-harmonic  point  force  given  by 


Re[F0exp(-iw0t) ] 


I*- 

v 


at  the  origin.  The  half -space  z>0  is  occupied  by  an  inviscid 
compressible  fluid  moving  with  constant  velocity  U“Mc  in  the 
x-direction,  where  M  is  a  dimensionless  Mach  number.  The  region 
z<0  is  in  a  vacuum.  The  Cartesian  co-ordinates  x,  y  and  z  are 
assumed  to  be  stationary  (i.e.  not  moving  with  the  fluid).  Figure 
1(a)  shows  the  geometry  of  interest. 

Let  w  ( x ,  y ,  t )  ,  p  ( x ,  y ,  z ,  t )  and  <p  ( x ,  y ,  z ,  t )  be  the 
instantaneous  values  of  plate  displacement,  acoustic  pressure  and 
acoustic  velocity  potential,  respectively.  Then  the  plate 
equation  of  motion  is 


D[ { ( 3  2/3x2 )  +  (32/3y2) }2  +  psh(a2/3t2)]w 
-  Re[F0exp(-iw0t) ] -0(x) -0(y)  -  p(x,y,0,t) 


(1) 


and  the  linearized  convected  form  of  the  acoustic  wave-equation 
is 


v2*>  -  (1/c2)  [  (3/at)+U(3/3x)  ]  2*> 


(2) 


The  boundary  condition  at  the  plate  surface,  reflecting 
continuity  of  displacement,  is 


( 3<p/3z )  (x  ,y ,  0 ,  t )  -  -p((3/3t)  +  U(3/3x)]w  (3) 


and  the  acoustic  pressure  p  can  be  expressed  in  terms  of  the 
potential  <p  as 


f.'.j 


p  -  -p(  (3/3t)  +  U(  3/3x)  ]*> 


(4) 


By  defining  complex 
<D(x,y,z)  of  w,  p  and 


amplitudes  W(x,y) 
<0  such  that 


P ( x , y , z )  and 


w(x,y, t) 


Re(W(x,y) • exp(-iw0t) ] 


(with  similar  expressions  for  p 
and  0  as  Fourier  transforms, 


and  <p) ,  and  expressing 


oo 

W(x,y)  -  (1/4tt2)  JJ  W(a,/9)exp(  iax+i/?y)  dad/9 


P(x,y,z) 

-  (1/47T2)  [ 

as 

P(a,/9,z) 

0(x,y,z) 

J 

oo 

i(a,/?,  z) 

exp(  iax+i/9y)  dad/? 


the  solution  for  these  variables  can  be  shown  to  be 


W(x,y)  -  (F 


Jl  — 


P(x,y,  z)  -  (-ipo/F  /4tt2) 

oo 


oo 

Tii 

»  r 


-aM/k)  exp(  iax+i/9y+j-yz) 


Z(wQ  ,a,/9) 


dad/9 


0(x,y,z)  =  (_waF0/47r2) 


(1-aM/k)  exp(  iax+ifly+j-yz ) 


Z(wo  ,a,/9) 


dad/9 


where 

Z(wo,a,£)  -  D(a2  +  B2)2  -  Pghw2  -  ipw2 ( 1-aM/k) 2 /y 
is  the  dispersion  relation  and 

y  m  V(k2(l-aM/k)2  -  a2  -  B2),  a2+/?2<k2(  1-aM/k)  2 

>  -  iV{a2  +  B2  -  k2  ( 1-aM/k)  2  } ,  a2+/?2  >k2  ( 1-aM/k)  2 


to  satisfy  the  radiation  condition  for  outgoing  waves. 


W,  P 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(ID 

(12) 
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The  acoustic  velocity  amplitudes,  Vx,  Vy  and  Vz  (or 
30/3x,  30/ 3y  and  30/3z)  ars  given  by 


vx  -  t-iUoP0 


/«-*>  jJ  9xp<^:^rz>  — * 


(’.3) 


Vy  "  (-iWoF0 


/4"!)  JJ  mLfiMggn.  dadfl  (14, 


V  =  (-iw 
z  o 


'./to1)  Jja-^)  d°d* 


(lb) 


3.  LINE-EXCITED  INFINITE  PLATE 


The  geometry  of  interest  is  shown  as  Figure  1(b).  The  plate 
displacement  W,  the  acoustic  pressure  P,  the  acoustic  potential 
0  and  the  particle  velocities  Vx  and  Vz  are  derived  via  the 
one-dimensional  Fourier  transform  as 


W(x)  -  (P,/2»)  T  do, 

-i  0 


(16) 


QO 


o  a' 


Vx  “  (-iwoPo/27r) 


y 

Z(wq , a) 

(1-asM/k) 

exp( iax+iyz) 

y 

Z(wQ,a) 

a( 1-aM/k) 

exp( iax+iyz ) 

(17) 


(18) 


Z(u  ,a) 


■y 


(19) 


(20) 


vz  -  (-i..F0/a„  Ja-OMA)  d. 

—to  ® 

where 


2(wo 

,a)  -  Da4  -  pghu2  - 

ipw2 (1-aM/k) 2 /y 

(21) 

y  = 

V{fc2 ( 1-aM/k) 2-a2 } , 

lal <kl 1-aM/k 1 

(22) 

y  - 

iV(a2-k2(l-aM/kfi, 

lal >kl 1-aM/k 1 

(23) 

4.  DEFINITION  OF  THE  ACOUSTIC  INTENSITY  VECTOR 


is  defined  as  the  time- 
intensity  or  acoustic 


(24) 


in  a  moving  acoustic 
medium  has  been  a  subject  for  debate  in  the  literature  [9-12]. 

Two  alternative  definitions  of  N  have  been  proposed,  both  of 
which  satisfy  a  conservation  equation,  viz. 

6E/at  +  div  N  -  0  (25) 

where  E  is  the  instantaneous  acoustic  energy  density.  The  first 
definition  (referred  to  by  using  subscript  1),  according  to 
Morfey  [10],  gives  the  following  expressions  for  Ex  and  the 
components  NXI,  Nyt  and  NZ1  of  Nj^ 


The  acoustic  intensity  vector  I 
average  of  the  instantaneous  acoustic 

energy  flux  vector  N. 


I  -  lim  t(l/T)J  N  dt] 

T-«o  -i 


The  problem  of  defining  the  vector  N 


E 


I 


fp((W>)2 


N. 


xi 


+  (1/C2)  (*>t+Uox)2  - 
-  -pv?t{  (l-U2/c2)<Px- 

Nyi  «  -POt<Py 


(2U/c2)^x(«Pt+U^x)} 

(U/c2)<Pt) 


(26a) 

(26b) 

(26c) 

(26d) 
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The  second  definition  (using  subscript  2),  according  to  Ryshov  & 
Shefter  [9],  is 

E2  =  \p{(v*)2  +  (1/c2)  (*>t+ *px)2}  (27a 

NX2  “  -P(*t+  U'px)'f>x  +  5PU{(V^)2  +  (1/c2)  («Pt+U^x)2}  (27b 

Ny2  *  -p(-Pt+U-Px)<Py  (27c 

NZ2  =  -P(‘Pt+U«Px)<Pz  (27d 

In  the  above  expressions,  subscripts  t,  x,  y  and  z  of  <p 
denote  partial  derivatives  with  respect  to  time  and  the  co¬ 
ordinate  directions  respectively. 

The  quantities  and  N ^  (i=l,2)  can  be  expressed  in 

terms  of  tne  instantaneous  acoustic  pressure  p  and  velocity 
components  vx,  Vy  and  v2  as 


Ei  =  p2/ ( 2pc2 )  +  (U/c2)pvx  +  ^p(v2+vy+v2) 

(28a 

Nxi  =  Pvx(l+M2)  +  (M/pc)p2  +  pUv2 

(28b 

Nyi  =  (P+PUvx)vy 

(28c 

Nzi  =  (P+PUvx)vz 

(28d 

E2  =  p2/ ( 2pc2 )  +  ^p(v2+vy+v2) 

(29a 

NX2  "  pvx  +  (M/2P°)P2  +  ;pU(v2+v2+v2) 

(29b 

N  =  pv 

y  2  v  y 

(29c 

N  =  pv 

Z  2  *  Z 

( 29d 

It  should  be  noted  that  the  expressions  for  the  components  of 
— ♦ 

Nx  and  N2  reduce  to  well-known  expressions  in  the  absence  of 
flow  (U-M-O). 


>-*•-  ft. 


i 


The  components  of  the  t*  s-averaged  intensity  vector  I 

according  to  formulations  1  and  2  are,  in  terms  of  the  complex 

amplitudes  P,  V  ,  V  and  V  , 

x  y  z 

lxi  "  ^Re[{P(l+M2)+pUVx}*Vx  +  (M/pc) P*P]  (3 

Iyi  "  (3 

!2i  =  ^Re[{P+pUVx}*Vz]  (3 

=  tRe[P*Vx  +  (M/2pc)P*  P  +  fpU(VxVx+V*Vy+V*Vz)]  (3 


I  »  -Re[ P*V  1 
y  2  2  1  y 


I  -  -Re[ P*V  ] 
22  2  1  2  J 


(31c 


In  the  case  of  a  line-excited  plate  without  flow,  these 
expressions  reduce  to  those  used  by  Spicer  [8]. 


The  far-field  sound  pressure  can  be  evaluated  for  both 
point-  and  line-force  excitation  using  the  method  of  stationary 
phase,  which  is  described  in  some  detail,  for  the  case  of  point 
force  excitation,  in  Appendix  A. 


The  far-field  pressure  in  the  case  of  point-force  excitation 
is  best  evaluated  using  a  non-standard  spherical  co-ordinate 
system  identical  to  that  used  by  Leppington  &  Levine  120]  to 
calculate  the  far-field  pressure  of  a  circular  piston  set  in  an 
infinite  rigid  baffle.  The  co-ordinate  system  is  defined  by 


x  -  R*cos0,  y  -  R-sinP'Siny,  z  *  R-sinP-cos? 


and  illustrated  in  Figure  2(a).  The  stationary  phase  point 


(<*„,£„)  is  given  by 


a0  -  [k/  ( 1-M2 )  ]  (-M  +  cose  /  V{h(M,0)}] 


0O  *  k’Sine-sinir  /  V{h(M,0)} 


where 


h(M,0)  -  1  -  M2s in2e 


The  resulting  expression  for  the  far-field  pressure  is 


P(R,e,*) 


P  [=AS](i-a  MA)2fain?i£§M 

o  2ffR  ' _ a  '  '  h(M,0) 


exp{ iR*g(k,M, 9) } 


/  -i  _  /i-x  2  f  ik*  sine  •  cosiri  , ,  ,  2  .  „2 „  2  „  4 , 
(l-«tM/k)  -  [„.-v(h(M,9)/jtl-(»0^0)  /kp) 


where 


g(k,M, 0)  -  [k/(l-M2)] [-M-cose  +  V( 1-M2s in20) ] 


The  expression  for  g(k,M,e)  is  identical  to  that  obtained  by 
Leppington  &  Levine  for  the  piston  problem.  It  can  be  shown  that 
equation  (36)  reduces  to  the  well-known  classical  result  in  the 
no-flow  case.  Equations  (32)  to  (37)  are  expressed  in  terms  of 
the  traditional  spherical  co-ordinate  system  (Figure  2(b))  in 
Appendix  B. 


A  cylindrical  co-ordinate  system  is  used  defined  by 


x  -  R’cos®,  y  “R’SinS 


and  illustrated  in  Figure  2(c).  The  stationary  phase  point  is 
given  by 

a0  -  [ k/ ( 1-M2 ) ] [ -M  +  cose  /  V{h (M, 0 ) } 1  ( 


V  '-' .  V  -*  l1' 


where  h(M,0)  is  given  by  equation  (35).  Equation  (39)  gives  an 
expression  for  a0  which  is  identical  to  that  in  the  case  of 
point-force  excitation,  equation  (33).  The  far-field  pressure  is 
given  by 


P(R,©) 


F 

0 

f-ik] 
!  2ttrJ 

1/2  ( 1-atgM/k)  2 

sin© 

h(M,0)3/4 

exp{ iR*g (k,M, 0) } 

(l-aflM/k)2  - 

[  ik-sin© 

[/i  •  V(h  (M,  © )  }  J 

(40) 


where  g(k,M,0)  is  given  by  equation  (37). 

Inspection  of  equations  (36)  and  (40)  reveals  that  the 
expressions  for  the  far-field  pressure  can  varnish  if  a0=k/M,  but 
this  is  only  true  if  M*l.  For  the  range  of  Mach  numbers 
encountered  in  water-loaded  plates  the  theory  indicates  that  the 
flow  has  only  a  very  small  effect  and  this  is  borne  out  in  the 
numerical  calculations  described  in  Section  7.3. 


6.  DISCUSSION  OF  STABILITY 

(a)  Ste.ady-3t^tg_ABPJ.P.ftth 


It  is  instructive  at  this  stage  to  apply  the  results  of 
Benjamin  (13,14]  and  Landahl  [15]  to  the  stability  of  a  line- 
excited  plate.  Their  results  are  obtained  for  the  general  problem 
of  inviscid  fluid  flow  over  a  flexible  solid  plane  boundary. 
Provision  is  made  for  including  damping  in  the  solid  and  the 
effect  of  a  thin  boundary  layer  in  the  fluid.  In  the  absence  of  a 
boundary  layer  and  for  an  undamped  plate,  the  equation  for  the 
phase  speed  can  be  expressed  in  the  present  notation  as 

(U-Cp) 2  -  (Pgha/p) (c2-c2)  (41) 

where  Cp  is  the  phase  speed  of  free  waves  in  the  coupled  plate- 
fluid  system  and  c0  the  speed  of  in  vacuo  free  waves  in  the 
plate,  given  by 


c o  -  a-V(D/psh) 


Wr iting 


C  -  Pgha/p 


l6- 


equation  (41)  can  be  expressed  as  a  quadratic  in  c0  which  can 
be  solved  to  give  c0  as  a  function  of  C  and  c0  (which  are 
both  functions  of  a),  viz. 


cp  -  [1/U+1)][U  +  VU(C+l)c2  -  <U2}]  (42) 

For  a  fixed  positive  value  of  a,  there  are  three  distinct 
possibilities,  which  Benjamin  and  Landahl  interpret  as  follows: 

2  2 

(i)  U  <Ccfl,  which  gives  one  positive  and  one  negative  root.  The 

negative  root  corresponds  to  a  wave  propagating  in  the  opposite 
direction  to  the  flow.  Benjamin  [14]  showed  that  this  wave  decays 
if  damping  is  present  in  the  plate.  The  wave  with  positive  Cp 
also  propagates  without  damping  in  the  plate,  but  it  may  be 
destabilized  in  the  presence  of  a  boundary  layer.  This  wave 
corresponds  to  a  class  B  instability,  as  mentioned  in  the 
Introduction. 

(ii)  <c2<U2< (C+l)c2 ,  which  gives  two  positive  roots,  one  less 

than  CJ  and  the  other  greater  than  U.  The  faster  wave  has 
properties  similar  to  the  wave  with  positive  cp  in  case  (i).  In 
a  uniform  flow,  the  slower  wave  is  destabilized  by  damping  in  the 
plate,  corresponding  to  a  class  A  instability.  The  transition 
between  cases  (i)  and  (ii),  i.e. 

U2  =  £c2  or  c  «  0 
o  p 

cam  be  expressed  as 

pU2  «  Da3 


This  corresponds  to  the  condition  known  as  "static  divergence"  in 
studies  of  panel  flutter.  The  hydrodynamic  pressure  is  just 
sufficient  to  maintain  a  stationary  wave  against  the  resilient 
restoring  forces.  For  an  infinite  plate,  this  condition  occurs  at 
zero  frequency  giving  an  infinite  wavelength,  as  mentioned  by 
Dowell  [21].  Dowell  pointed  out  that  no  physical  significance  can 
be  attached  to  this  result  for  an  infinite  plate,  but  for  a  plate 
of  finite  width  or  an  infinite  plate  resting  on  an  elastic 
foundation,  both  the  critical  wavelength  and  flutter  velocity 
become  finite. 


U2>U+l)c2 


which  gives  two  complex  values  of  c  ,  one  of 


(iii) 

which  has  a  positive  imaginary  part  and  thus  represents  a 
propagating  wave,  the  amplitude  of  which  grows  exponentially  with 
time.  This  is  a  Class  C  instability. 


The  above  approach  to  stability  is  compared  with  that  using 
the  causal  approach  (see  Section  6(b)  below)  in  Section  8. 
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In  Sections  2  to  5  it  is  assumed  that  the  excitation 
F0exp(-i«0t)  exists  for  all  time  t.  In  practice,  however,  the 
excitation  may  be  switched  on  at  some  instant,  say  t-0,  and 
steadily  maintained  thereafter.  If  this  is  the  case,  then  there 
are  three  possibilities  for  the  response  at  any  point  of  the 
system  after  a  long  time  has  elapsed.  The  amplitude  may  either 
(i)  settle  down  to  the  steady-state  value,  (ii)  settle  down  to  a 
different  steady-state  value  but  increase  exponentially  with 
distance,  or  (iii)  grow  exponentially  with  time  at  every  point  of 
the  system.  Case  (ii)  is  termed  a  convective  instability  and  case 
(iii)  an  absolute  instability  (Melcher  [18]). 

This  form  of  excitation  can  be  represented  mathematically  by 
a  Laplace-Four ier  transform.  To  simplify  the  analysis,  the 
specific  example  of  the  displacement  of  a  line-excited  plate  is 
discussed.  However,  the  technique  described  here  can  be  applied 
to  the  plate  displacement,  the  acoustic  pressure  and  the  acoustic 
particle  velocities  of  both  point-  and  line-excited  plates. 

Suppose  that  the  plate  displacement  w(x,t)  can  be 
expressed  in  terms  of  a  complex  time-dependent  amplitude 


w(x,t)  -  Re[W(x, t) ] 


Then  the  Laplace-Four ier  representation  of  W  is 


«»+ia 


W(x,  t)  =  (1/2 it)  J  W(x,«j)  exp(-iut) 


doj 


(43) 


-oo+icr 


W(x,cj)  -  (1/2tt)  W(a,u)  exp(iax)  da 


(44) 


where 


W(a,u)  -  iF 0/{  (w-cjq  )  Z  ( a,  w)  }  (45) 

and  a  is  a  positive  constant  such  that  the  contour  Im(u)*a 

lies  above  any  singularities  of  W(x,u).  This  ensures  that 
causality  is  satisfied,  i.e.  that  the  response  of  the  system  is 
zero  for  t<0.  In  equation  (45),  Z(a,w)  is  given  by  equation 
(21)  with  u0  replaced  by  w. 


Provided  that  the  radiation  condition  is  satisfied  along  the 


real  axis,  the  integrand  of  equation  (44)  can  be  analytically 
continued  into  the  complex  a-plane.  However,  the  integrand  has 
branch  points  at 

a  *  a~  -  -k/(l-M),  a  »  a+  =  k/(l+M) 


corresponding  to  7=0.  For  a>0,  these  branch  points  are  off  the 
real  a-axis  and  in  opposite  half -planes.  The  associated  branch 
cut  must  not  cross  the  real  a-axis  and  must  therefore  go  to 
infinity,  the  cut  from  a+  going  to  +«  and  that  from  a~ 
going  to  -<*».  Apart  from  this  restriction,  the  choice  of  branch 
cuts  is  arbitrary.  For  simplicity,  the  branch  cuts  used  in  all 
the  numerical  work  are  defined  by  Im[7]=0,  as  a  result  of  which 


7  =  V{k2(l-aM/*)2  -  <*2},  Im[7]>0  (46) 


everywhere,  except  on  the  branch  cuts. 

The  Fourier  integral  (44)  can  be  evaluated  by  suitably 
closing  the  contour  and  applying  Cauchy's  theorem.  Figure  3 
illustrates  the  case  when  x>0,  requiring  the  contour  to  be 
closed  in  the  upper  half -plane.  For  x<0,  the  contour  is  closed 
in  the  lower  half -plane.  This  convention  ensures  that,  as  Im[a] 
tends  to  infinity,  the  term  exp(iax)  in  the  integrand  remains 
finite.  In  Figure  3,  let  a-j  +  (j-1,3)  represent  all  the 
singularities  of  the  integrand  in  the  upper  half -plane  and  a-j- 
(j=l,3)  those  in  the  lower  half -plane.  If  a  is  large  enough, 
then  all  the  poles  are  well  clear  of  the  real  axis.  The  contours 
Cx  and  C2  are  circular  arcs  of  large  radius  and,  provided  that 
all  singularities  in  the  upper  half -plane  lie  inside  the  contour, 
the  contributions  from  C2  and  C2  tend  to  zero  as  the  radius 
tends  to  infinity.  The  integrals  along  Blf  B2  and  B3  are 
branch  cut  integrals  which  are  discussed  later  in  the  Section. 

The  integral  (44)  can  be  expressed  using  Cauchy's  theorem  as  the 
sum  of  residues  and  branch  cut  integrals,  in  the  form 


-v  i  i 

W(x,td)  «  ±  E  -F  exp(  ia .  x) /{ (w-u  )Z'(a.,w)}  +  branch  cut  integrals 
j-l  °  J  °  1  (47) 


where  the  +  sign  corresponds  to  positive  values  of  x  and  the 
-  sign  to  negative  values  of  x. 

The  plate  response  is  given  by  using  equation  (43)  to 

integrate  W(x,w)  over  the  Laplace  contour.  For  t<0,  the  contour 
must  be  closed  as  in  Figure  4(a)  to  be  consistent  with  causality. 


For  t>0,  any  numerical  scheme  to  evaluate  the  integral  (43) 
would  eventually  become  unstable  at  large  t  because  the  term 
exp(-iart)  in  the  integrand  tends  to  infinity  as  t-«».  For  this 
reason  and  also  to  be  able  to  provide  a  physical  interpretation 
of  the  response  when  t  is  large  and  positive,  it  is  necessary 
to  move  the  Laplace  contour  down  to  the  real  axis,  i.e.  to  let  a 

tend  to  zero.  If  there  are  no  singularities  of  W(x,cj),  apart 
from  that  at  cj*w0  ,  then  the  Laplace  contour  for  large  t  can  be 
moved  so  that  it  lies  just  below  the  real  axis,  indented  around 
(j*u0  and  closed  in  the  lower  half-plane,  as  in  Figure  4(b).  The 
plate  response  can  now  be  expressed  as  the  sum  of  modes  and 
branch  cut  integrals 


W(x,t)  =  £  W*  exp(-a*x)exp(-iw  t)  +  branch  cut  integrals  (48) 

j-l  3  3  ° 

where 

W*  =  ±{-Fo/2' (a*,wQ)} 

In  order  to  ascertain  whether  or  not  any  physical 
s ignif icance  can  be  attached  to  each  mode  of  the  plate  response, 
it  is  necessary  to  use  the  root  locus  technique  described  by 
Melcher  [18].  First  consider  when  there  are  no  singularities  of 

W(x,cj)  in  the  upper -half  w-plane,  and  ignore,  for  the  moment, 
the  contribution  of  the  branch  cut  integrals.  For  a  given  cj0  , 
the  loci  of  the  roots  aj*  (which  are  functions  of  u)  are 
plotted  for  values  of  cj=cjQ  +  io  as  a  varies  from  +«>  to  0.  As 
the  Laplace  contour  is  moved  down  to  the  real  u-axis,  the  poles 
enclosed  by  the  Fourier  contour  also  move  and  it  may  be  necessary 
to  deform  the  Fourier  contour  so  that  it  always  includes  the  same 

poles  of  W(a,u).  There  are  six  basic  possibilities,  illustrated 
in  Figure  5(a).  The  shaded  circles  in  this  Figure  denote  points 
where  the  loci  terminate,  corresponding  to  w«u0. 

Suppose  that  a  locus  comes  from  infinity  in  the  upper  half 
a-plane  to  a  value  above  the  real  axis,  say  ar  +  iOi  with  a^X) 
(case  (i)  in  Figure  5(a)).  The  response  of  this  mode  is  nonzero 
only  for  x>0,  in  which  case 


W(x,t)  -  Wj^expt  (  iar-ai)x}exp(-iw0t) 


which  gives  a  decaying  or  "evanescent"  wave.  If  the  locus 
terminates  on  the  real  axis  (case  (ii)),  the  response  for  x>0 
becomes 


W(x,t)  =  W2exp(iarx)exp(-iu0t) 


which  represents  a  travelling,  non-decaying  wave.  Since  the  locus 
comes  from  above  the  real  axis,  the  Fourier  contour  must  be 
indented  below  the  pole.  The  wave  is  travelling  in  the  x- 
increasing  direction,  (i.e.  outwards  from  the  excitation)  if  ar 
is  positive  and  in  the  x-decr easing  direction  (i.e.  inwards  from 
infinity)  if  «r  is  negative.  If  the  locus  crosses  the  real  a- 
axis  (case  (iii)),  the  Fourier  contour  can  no  longer  be  taken 
along  the  real  axis  but  must  be  deformed  around  the  pole.  If  the 
pole  is  at  a*ar-ia^,  where  a^>0,  then  the  response  for  x>0  is 


W(x, t)  -  W3 exp{ ( iar+a 2 )x} exp(-iw0t) 


which  for  positive  ac  denotes  a  wave  propagating  in  the  x- 
incr easing  direction  and  increasing  exponentially  in  amplitude 
with  distance.  This  is  usually  termed  a  "convective  instability". 
Cases  ( iv)  to  (vi)  are  identical  to  cases  (i)  to  (iii) 
respectively,  except  that  the  responses  are  non-zero  only  for 
x<0 . 


Suppose  there  are  zeros  of  Z'(aj  ,<j)  in  the  w-plane,  where 
aj  is  one  of  the  the  roots  of  the  dispersion  relation  Z(a,u)*0. 
Then  these  must  be  identified  as  possible  singularities  of 

W(x,gj).  Melcher  [18]  shows  that  a  singularity  of  W(x,u)  in  the 
w-plane  is  the  result  of  the  coalescence  of  a  pair  of  poles  in 
the  a-plane.  He  also  shows  that  there  are  four  distinct 
possibilities  here,  illustrated  as  Figure  5(b).  Suppose  that  both 
loci  come  from  either  above  or  below  the  real  a- ax is  and 
coalesce  before  reaching  the  axis  (case  (i))  or  on  the  axis 
(cases  (ii)  and  (iii)).  Melcher  points  out  that  the  terms  of 

W(x,w)  corresponding  to  each  pole  contain  a  singularity,  but 
when  they  are  combined  the  singularities  are  not  present.  Hence 
the  only  singularity  in  the  u-plane  is  at  a>”w0  and  the 
Laplace  integral  can  be  evaluated  giving  equation  (48) .  The 
physical  signif icance  of  the  resulting  modes  then  depends  on 
where  the  loci  terminate,  as  shown  in  Figure  5(b).  In  case  (iv), 
a  locus  coming  from  above  the  a-axis  coalesces  with  one  coming 
from  below.  In  this  case,  Melcher  showed  that  the  Fourier  contour 
can  no  longer  be  distorted  to  include  the  same  poles  but  will 
always  include  one  or  the  other  when  closed  in  the  upper  or  lower 
half  a-plane.  In  the  <J-plane  there  is  a  branch  pole  in  the 
upper  half-plane,  at  cj-cj,  ,  and  when  the  Laplace  contour  is 
deformed  around  this,  there  is  always  a  contribution  exp(ust) 
to  the  response,  with  a>3>0.  Hence  the  response  grows 
exponentially  with  time  at  a  fixed  location  x.  This  is  termed  an 
"absolute  instability",  which  occurs  whatever  the  value  of  u0 . 


It  is  now  necessary  to  discuss  the  branch  cut  integrals.  In 


the  absence  of  flow,  Strawderman  et  al.  [6]  show  that  the 
contribution  arising  from  the  circular  contour  B2  around  the 
branch  point  is  zero,  and  that,  provided  the  branch  cuts  are  well 
away  from  any  singularities  of  the  integrand,  the  integrals  along 
and  B3  contribute  only  to  the  near  field  because  the 
integrands  are  bounded.  It  can  be  shown,  by  a  simple  linear 
mapping,  that  the  integrals  with  flow  can  be  converted  to  the 

same  form  as  those  without  flow,  and  hence  that  the  above  results 

also  apply  in  the  presence  of  flow.  Crighton  [7]  has  pointed  out, 
for  the  case  of  an  undamped  plate  without  flow,  that  it  is  only 
those  poles  that  exist  for  all  suitable  choice  of  branch  cuts 
which  have  any  physical  significance  at  large  x.  The  remaining 
poles  contribute  only  to  the  near  field  and  their  contributions 
must  always  be  combined  with  those  from  the  branch  cut  integrals. 

As  already  mentioned  in  this  Section,  if  a  is  large 
enough,  all  the  poles  of  the  integrand  of  equation  (44)  are  well 
away  from  the  real  axis.  A  consequence  of  this  is  that  the 
location  of  the  poles  depends  on  the  particular  choice  of  branch 

cuts  and  that  no  pole  exists  for  all  possible  choices  of  branch 

cuts.  It  also  means  that,  as  o  varies,  a  root  locus 
corresponding  to  the  movement  of  a  particular  pole  may  cross  a 
branch  cut  and  "disappear",  or  vice  versa.  At  o=0 ,  the  poles 
corresponding  to  evanescent  waves  also  depend  on  the  choice  of 
branch  cuts.  The  contribution  to  the  response  arising  from  these 
poles  can  therefore  have  no  physical  significance  in  isolation, 
but  must  be  combined  with  the  branch  cut  integrals.  A  free-wave 
pole  (i.e.  one  on  the  real  axis  for  a»0)  always  exists  for  all 
possible  choices  of  branch  cut  and  therefore  has  physical 
significance.  A  pole  corresponding  to  a  convective  instability 
crosses  the  real  axis  at  some  positive  value  of  a  and  always 
exists  in  a  neighbourhood  of  the  real  axis.  Furthermore,  if 
coalescing  poles  in  the  a-plane  corresponding  to  a  singularity 
in  the  upper-half  w-plane  exist,  they  do  so  for  all  possible 
choices  of  branch  cut. 

^  To  summarize,  provided  that  there  are  no  singularities  of 
W(x,u)  in  the  w-plane,  the  root  locus  technique  cam  be  used  to 
determine  the  physical  significance  of  a  wave  having  complex 
values  of  a  for  real  ui.  The  modes  of  the  response  cam  be  either 
non-decaying  travelling  waves  or  convectively  amplifying  waves. 

If  a  singularity  exists  anywhere  ip  the  upper  half  w-plane  as  a 
result  of  coalescing  poles  whose  loci  in  the  a-plane  come  from 
opposite  half -planes,  then  the  system  has  an  absolute  instability 
and  the  response  grows  with  time  at  a  fixed  location.  In 
practice,  this  is  usually  damped  by  some  nonlinear  mechamism  or  a 
factor  not  taken  into  account  in  the  current  theory.  A  necessary 
condition  for  instability  (absolute  or  convective)  is  that  a 
locus  in  the  a-plane  crosses  the  real  axis.  This  implies  that 
for  some  real  wavenumber  a,  the  dispersion  relation  must  have  a 
root  w  with  a  positive  imaginary  part,  the  criterion  used  by 
Benjamin  [13,14],  Landahl  [15]  and  Miles  [16]  to  determine 
instability. 


-22- 


a*i;i 


fvf 


iff 


All  the  numerical  results  presented  in  this  section  are  for 
an  infinite  steel  plate  of  thickness  5  cm,  with  E=19*5xl010  N/m2, 
v=*0*29  and  pg=7700  kg/m3.  The  plate  is  immersed  in  water  of 
density  1000  kg/m3  and  sound  speed  1500  m/s.  The  quantity  F0  is 
IN  in  the  case  of  point-force  excitation  and  1  N/m  in  the  case  of 
line-force  excitation.  Plate  damping  when  used  is  represented  by 
a  complex  Young's  modulus  E=E(l-i7j),  where  7}  is  a  hysteretic 
loss  factor.  Two  flow  speeds  are  considered,  U=7 • 5  and  22*5  m/s 
(15  and  45  knots),  equivalent  to  Mach  numbers  of  0*005  and  0*015 
respectively.  Most  of  the  numerical  results  are  limited  to  the 
case  of  line-excited  plates,  but  the  far-field  pressure 
computations  are  for  both  point-  and  line-excited  plates. 


For  an  undamped  line-excited  plate,  the  real  roots  of  the 
dispersion  relation  Z(cj0,a)-0  (where  Z(u0,a)  is  defined  by 
equation  (21))  determine  the  wave  numbers  a  corresponding  to 
free  waves  propagating  along  the  plate  at  a  given  frequency  u0 . 
Positive  values  of  a  correspond  to  waves  propagating  in  the  x- 
increasing  direction  and  negative  values  to  waves  propagating  in 
the  x-decreasing  direction. 

Figure  6  shows  a  plot  of  wavenumber  a  against  frequency  f 
(~u0/2tt)  for  a  line-excited  plate  with  a  7*5  m/s  flow,  which  is 
compared  with  the  corresponding  plot  in  the  absence  of  flow.  At 
low  frequencies  the  flow  has  a  significant  effect  on  wavenumber 
at  a  given  frequency.  Further  results  not  presented  here  indicate 
that  with  increasing  frequency  the  values  of  a  with  flow 
approach  those  without  flow.  At  frequencies  above  coincidence 
(approximately  5  kHz),  values  of  a  with  and  without  flow  are 
almost  identical.  In  the  case  with  flow,  at  frequencies  between 
fL  (approximately  0*109  Hz)  and  fy  (approximately  0*121  Hz),  the 
number  of  positive  roots  of  the  dispersion  relation  changes  from 
one  to  three.  At  frequencies  fL  and  fg,  two  of  these  positive 
roots  coalesce.  This  effect  is  discussed  in  Section  7(d)  below. 

At  a  frequency  f q  of  approximately  0*112  Hz,  the  curve  with 
flow  for  positive  values  of  a  has  a  local  minimum  value. 

Figure  7  shows  plots  of  phase  velocity  wQ/a  and  group 
velocity  dua/da  against  frequency  at  a  flow  speed  of  7*5  m/s.  At 
frequencies  below  fL  (0*109  Hz),  the  positive  value  of  phase 
velocity  is  less  than  the  flow  speed  while,  at  frequencies  above 
fjj  (0*121  Hz),  it  is  greater  than  the  flow  speed.  At  frequencies 
above  coincidence,  further  numerical  results  not  presented  here 
indicate  that  the  phase  speed  of  the  upstream  wave  tends  to  c-U 
and  that  of  the  downstream  wave  to  c-KJ.  At  the  frequency  f^ 
(0*112  Hz),  corresponding  to  the  local  minimum  in  the  wavenumber- 


frequency  curve  (Figure  6),  the  group  velocity  has  a  singularity. 
At  low  frequencies,  there  are  two  negative  values  of  group 
velocity  which,  as  the  frequency  tends  to  zero,  both  tend  to  the 
same  limiting  value,  approximately  (-3/2)0.  This  limiting  value 
can  also  be  obtained  from  the  dispersion  relation  by  expressing 
it  as  a  function  of  a  and  u0  fund  talking  the  limit  as  w0-Q 
of 


-(dZ/da)/(dZ/du0) 


The  result  is  (-3/2)U  multiplied  by  a  negligibly  small  Mach 
number  correction  arising  from  the  compressibility  of  the  fluid. 
If  the  plate  is  substituted  by  a  membrane,  with  Da4  in  the 
dispersion  relation  replaced  by  Ta2 ,  where  T  is  the  tension  in 
the  membrane,  then  the  limiting  group  velocity  at  zero  frequency 
becomes  (-1/2)U  with  the  same  Mach  number  correction.  This 
effect  is  discussed  below  in  Section  8.  This  corresponds  to 
static  divergence  as  mentioned  in  Section  6(a). 

Figure  8  shows  a  plot  of  wavenumber  against  frequency  at  a 
flow  speed  of  22*5  m/s,  which  is  also  compared  with  the 
corresponding  no-flow  plot.  The  dispersion  relation  only  has  one 
real  positive  root  a  over  the  whole  frequency  range,  which 
implies  that  there  are  no  real  values  of  a  and  for  which 

Z'(a,w0)=0,  and  hence  no  coalescing  poles  on  the  real  a-axis. 
Figure  9  shows  plots  of  phase  and  group  velocities  at  the  same 
flow  speed.  There  is  a  singularity  in  the  group  velocity  and  the 
limiting  value  of  group  velocity  at  zero  frequency  is  (-3/2)LJ, 
the  same  result  as  with  a  flow  speed  of  7*5  m/s. 

It  is  possible  to  find  all  the  complex  roots  a  of  the 
dispersion  relation  Z(a,ui)=0  corresponding  to  real  values  of 
ui.  However,  the  number  and  location  of  the  complex  roots  depends 
on  the  choice  of  branch  cuts,  and  therefore  (according  to 
Cr ighton  [7]),  no  physical  significance  can  be  attached  to  any  of 
the  complex  roots.  Another  approach  is  to  find  the  real  and 
complex  roots  u  corresponding  to  real  values  of  a,  since  this 
determines  free  waves  of  the  form  exp( iax-iwt)  with  a  real.  If 
lm(u]>0,  this  indicates  a  wave  the  amplitude  of  which  grows  with 
time.  Without  loss  of  generality,  Re[w]  can  be  restricted  to 
positive  values,  which  gives,  as  before,  positive  values  of  a 
corresponding  to  waves  propagating  in  the  x-increasing  direc¬ 
tion,  and  negative  values  to  waves  propagating  in  the  x- 
decreasing  direction. 

Figure  10  shows  wavenumber  plotted  against  the  real  and 
imaginary  parts  of  frequency,  at  a  flow  speed  of  7-5  m/s.  At 
values  of  a  greater  than  about  0 • 1  or  less  than  about  -0*29, 
the  dispersion  relation  gives  only  real  values  of  w 
corresponding  to  propagating  waves  identical  to  those  plotted  in 
Figure  8.  At  small  positive  values  of  a  up  to  about  0*1,  the 


dispersion  relation  gives  a  pair  of  complex  conjugate  roots  u. 
The  root  with  positive  imaginary  part  corresponds  to  an 
amplifying  wave  propagating  in  the  x-increasing  direction.  Figure 
11  is  similar  to  Figure  10,  but  is  at  a  flow  speed  of  22*5  m/s. 
This  also  gives  an  amplifying  wave  propagating  in  the  x- 
increasing  direction  for  positive  values  a  up  to  about  0*24.  It 
is  possible  to  deduce  from  Figures  10  and  11,  using  arguments  put 
forward  in  Section  6,  that  an  instability  of  some  sort 
(convective  or  absolute)  exists  at  both  flow  speeds.  However,  it 
is  not  possible  to  make  any  further  deduction  as  to  the  exact 
nature  of  the  instability  by  using  the  information  supplied  by 
Figures  6  to  11  alone. 


Some  numerical  comparisons  have  been  made  of  far-field 
pressure  with  and  without  flow  for  damped  point-  and  line-excited 
plates  with  a  loss  factor  7j«0*02  over  a  frequency  range  of  0*01 
Hz  to  10  kHz.  Corresponding  results  with  and  without  flow  are 
identical  to  within  0*1  dB  in  amplitude  and  2°  in  phase  at  all 
frequencies  up  to  1  kHz  (i.e  at  frequencies  below  coincidence). 
Tables  I  and  II  illustrate  the  differences  with  and  without  flow 
at  10  kHz  (i.e.  above  coincidence)  for  point-  and  line-force 
excitation  respectively.  The  far-field  sound  pressures  are 
calculated  at  a  nominal  radius  of  1  m  for  various  values  of  0 
(with  ir- 0  in  the  case  of  point-force  excitation)  and  the 
amplitudes  are  expressed  in  dB  re  1  fiPa.  The  comparisons  show 
that  the  main  effect  of  the  flow  is  to  slightly  change  the 
direction  of  the  coincidence  lobes,  while  not  changing  the  peak 
amplitude.  In  the  case  of  point-force  excitation,  results  with 
?h*90o  not  presented  here  show  much  less  of  a  flow  effect.  This 
indicates  that  the  predominant  flow  effect  is  felt  directly 
upstream  and  downstream  of  the  force  rather  than  in  any  other 
direction. 


Intensity  vectors  have  been  calculated  numerically  for  a 
line-excited  damped  plate,  by  adapting  the  computer  programs 
written  by  Spicer  [8]  to  incorporate  the  effect  of  flow.  Both 
formulations  of  acoustic  intensity  (equations  (30)  and  (31))  have 
been  incorporated  into  the  new  programs.  Figures  12  to  15 
illustrate  the  variation  with  frequency  of  intensity  vectors 
along  the  plate  surface.  In  these  Figures,  the  length  of  each 
vector  is  proportional  to  the  square  root  of  its  magnitude.  Each 
row  of  intensity  vectors  was  computed  separately  and  each  vector 
in  a  particular  row  war  normalised  with  respect  to  the  maximum 
value  in  that  row.  Figure  12  shows  intensity  vectors  according  to 
formulation  1,  at  a  flow  speed  of  7*5  m/s.  At  frequencies  above 
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and  below  0*12  Hz,  the  free  wave  propagating  away  from  the  force 
in  the  x- increasing  direction  corresponds  with  the  positive 
wavenumber  and  that  in  the  x-decr easing  direction  with  the 
negative  wavenumber.  At  intermediate  frequencies,  the  predominant 
acoustic  energy  flow  is  aligned  with  the  fluid  flow.  Figure  13 
shows  results  using  formulation  2  at  some  of  the  selected 
frequencies.  There  are  slight  differences  between  the  two  sets  of 
results.  Downstream  of  the  excitation,  both  formulations  give 
similar  results,  but  upstream,  at  frequencies  below  0*12  Hz,  the 
two  sets  of  results  are  approximately  180°  out  of  phase.  Figure 
14  compares  results  with  and  without  flow.  The  flow  is  seen  to 
have  an  effect  only  at  very  low  frequencies.  In  Figure  15, 
results  using  formulation  1  at  a  flow  speed  of  22*5  m/s  show 
similar  effects  to  those  in  Figure  12.  These  results  are 
discussed  further  in  Section  8. 

According  to  equations  (8),  (13),  (14)  and  (15),  at  a 
frequency  at  which  (l-aM/k)  vanishes,  then  the  contribution  to 
the  acoustic  intensity  vector  from  the  corresponding  free  wave 
pole  a  also  vanishes.  This  frequency  corresponds  to  that  where 
the  phase  speed  equals  the  flow  speed  and  it  is  given  by 


w0=U2V(Psh/D) 


corresponding  to  approximately  0*12  Hz  and  1*06  Hz  at  flow  speeds 
of  7*5  and  22*5  m/s  respectively.  At  7*5  m/s  (Figures  12  and  13), 
the  acoustic  intensity  at  0*12  Hz  does  not  decay  rapidly 
downstream  of  the  excitation,  whereas  it  does  so  at  0*14  Hz.  At 
22*5  m/s  (Figure  15),  the  downstream  acoustic  intensity  decays 
rapidly  at  1*0  Hz.  At  both  flow  speeds,  the  rate  of  decay  falls 
off  as  the  frequency  is  increased. 


(d)  Numerical  Determination  of  Stability  Characteristics  Using 
Causal  Approach 

In  order  to  use  the  root  locus  technique  to  determine 
stability,  the  complex  roots  a  of  the  dispersion  relation  for  a 
real  or  complex  frequency  u  need  to  be  found.  The  numerical 
procedure  used  to  do  this  was  to  apply  Newton's  method  in  the 
complex  plane  starting  from  each  point  on  a  grid.  Numerical  tests 
for  convergence,  divergence  and  finding  the  same  root  more  than 
once  were  incorporated  into  the  method.  At  low  frequencies,  the 
compressibility  of  the  fluid  can  be  neglected  and  the  dispersion 
relation  reduces  to  a  polynomial,  in  which  case  standard 
numerical  techniques  are  available  for  extraction  of  the  roots. 
Comparison  of  both  methods  yielded  almost  identical  results  up  to 
at  least  10  Hz. 


Figure  16  shows  root  loci  diagrams  at  a  flow  speed  of  7*5 
m/s  for  an  undamped  line-excited  plate.  The  loci,  plotted  for 


various  values  of  u0,  represent  the  movement  of  the  roots  a  of 
the  dispersion  relation 


Z (w,a)  -  0  ,  (j  =  (J0  +  io 


as  a  varies  from  a  large  positive  value  to  zero.  On  Figure  16, 
the  circles  represent  points  where  a=0.  All  the  loci  labelled  1 
correspond  to  a  frequency  u0/2n  of  0*04  Hz.  The  loci  labelled 
la  and  lb  terminate  on  the  real  a-axis  and,  using  the  results  of 
Section  6(b),  represent  outgoing  and  incoming  waves  respectively 
in  the  region  x<0.  The  locus  labelled  lc  crosses  the  real  a- 
axis  and  terminates  in  the  lower  half-plane.  This  represents  an 
outgoing  wave  in  the  region  x>0,  the  amplitude  of  which 
increases  exponentially  with  distance  (i.e.  a  convective 
instability) .  The  locus  labelled  Id  terminates  in  the  upper  half- 
plane.  If  a  is  greater  than  a  certain  value,  then  this  locus 
crosses  the  branch  cut  (which  at  low  flow  speeds  and  relatively 
small  values  of  a  is  almost  coincident  with  the  positive 
imaginary  axis)  and  "disappears".  The  locus  marked  le  only  exists 
for  relatively  large  values  of  a.  As  the  loci  Id  and  le  do  not 
exist  in  some  neighbourhood  of  the  real  a-axis,  then  (as 
mentioned  in  Section  6(b) ) ,  they  are  of  no  physical  significance 
in  isolation,  but  the  contribution  from  these  poles  must  be 
combined  with  the  branch  cut  integrals. 

At  a  frequency  of  0-08  Hz,  the  loci  (labelled  2a-2e  in 
Figure  16)  show  the  same  behaviour  as  those  at  0*04  Hz.  At  0*12 
Hz,  there  are  four  loci  (labelled  3a-3d)  which  terminate  on  the 
real  axis.  Reference  to  Figure  6  indicates  that,  at  this 
frequency,  the  dispersion  relation  has  four  real  roots.  There  are 
two  propagating  waves,  one  incoming  and  one  outgoing,  in  x<0 
and  two  outgoing  waves  in  x>0.  At  0*16  Hz,  two  of  the  loci  (4a 
and  4c)  terminate  on  the  real  axis  and  hence  represent 
propagating  outgoing  waves,  one  in  x>0  and  one  in  x<0.  The 
other  two  loci  (4b  and  4d)  represent  evanescent  waves  and  are  of 
no  physical  significance  in  isolation.  The  loci  at  0-2  Hz  show 
the  same  behaviour  as  those  at  0-16  Hz. 

The  above  results  indicate  that  there  are  three  distinct 
frequency  bands  where  the  results  have  different  interpretations. 
The  boundaries  between  these  bands  are  the  frequencies  fL  and 
fu  (approximately  0-109  and  0-121  Hz  respectively)  mentioned  in 
Section  7(a)  and  Figure  6,  where  there  are  coalescing  poles  on 
the  real  a-axis.  Melcher  [18]  shows  that  neither  case  of 
coalescing  poles  corresponds  to  an  absolute  instability.  Further 
numerical  work  not  presented  here  indicates  that  there  are  no 
other  cases  of  coalescing  poles  in  the  a-plane.  Although  there 
are  examples  in  Figure  16  of  loci  at  the  same  value  of  w0  which 
intersect,  the  points  of  intersection  correspond  to  different 
values  of  o  and  so  they  are  not  examples  of  coalescing  poles. 
Hence  there  is  no  absolute  instability  and  the  interpretation  of 


the  results  in  the  two  previous  paragraphs  is  valid.  This 
interpretation  is  compared  with  that  using  the  Benjamin-Landahl 
approach  to  stability  in  Section  8. 

The  numerical  results  presented  in  Figure  16  were  repeated 
with  plate  damping.  The  plate  loss  factor  chosen  was  0 *2,  rather 
higher  than  encountered  in  practice  in  order  to  exaggerate  the 
effects  of  damping  for  clarity.  The  results  are  shown  in  Figure 
17.  The  loci  corresponding  to  propagating  waves  with  negative  a 
move  downwards,  i.e.  away  from  the  axis  and  hence  represent  waves 
propagating  away  from  the  excitation  in  the  region  x<0  and 
slowly  decaying  to  infinity.  At  frequencies  below  fL  (0*109 
Hz),  the  locus  corresponding  to  the  incoming  wave  from  infinity 
in  x<0  now  crosses  the  real  axis,  which  suggests  that  this  wave 
is  destabilized  by  plate  damping.  Further  numerical  work  is 
necessary  to  clarify  and  quantify  this  effect.  At  frequencies 
above  fu  (0*121  Hz),  both  the  loci  corresponding  to  propagating 
waves  move  away  from  the  real  axis  and  hence  both  waves  slowly 
decay  to  infinity.  In  the  intermediate  region  corresponding  to 
four  propagating  waves,  all  the  loci  move  away  from  the  axis. 

Figure  18  is  similar  to  Figure  16  but  is  at  a  flow  speed  of 
22*5  m/s.  At  a  frequency  of  between  0*8  and  0*85  Hz,  the  diagram 
shows  coalescing  poles  at  a  point  not  on  the  real  a-axis,  one 
coming  from  above  and  one  from  below.  As  mentioned  in  Section  6, 
this  is  indicative  of  an  absolute  instability  where  the  linear 
theory  predicts  that  the  response  at  every  point  of  the  system 
grows  exponentially  with  time.  Further  numerical  result.*,  not 
presented  here  suggest  that  plate  damping  does  not  remove  the 
singularity,  but  merely  changes  its  location  slightly. 


8.  EISCU5SI.QN.  QEJBESULKa 


Comparison  of  the  causal  approach  and  the  Benjamin-Landahl 
approach  to  stability  indicates  that  provided  the  causal  approach 
does  not  give  an  absolute  instability,  both  methods  give  the  same 
interpretation  of  the  various  free  waves.  The  comparison  is 
illustrated  in  Figure  19.  It  should  be  noted  that  where  a 
convectively  unstable  wave  exists  it  will  dominate  the  response 
in  the  region  x<0  or  x>0.  The  amplitude  of  the  disturbance 
(assumed  by  the  linear  theory  to  be  infinitesimal)  increases  with 
distance  until  it  can  no  longer  be  assumed  to  be  infinitesimal, 
when  linear  theory  breaks  down  and  the  disturbance  is  damped  by  a 
non-linear  mechanism. 

The  next  question  to  be  resolved  is  the  interpretation  of 
the  intensity  vector  plots  in  the  light  of  the  stability  results. 
The  intensity  vector  plots  were  obtained  by  integrating 
numerically  along  the  real  a-axis.  With  damping  the  free  wave 
poles  in  the  a-plane  move  off  the  real  axis,  and  integrating 
along  the  real  axis  implies,  according  to  Crighton  [7],  that  the 
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free  wave  poles  always  represent  outgoing  waves,  as  is  the  case 
for  a  line-excited  plate  without  flow.  In  order  to  satisfy 
causality  in  systems  with  a  convective  instability,  the  Fourier 
contour  (in  the  a-plane)  must  be  deformed  around  the  pole  giving 
rise  to  the  instability,  and  the  plate  response,  acoustic 
pressure  and  acoustic  particle  velocities  must  be  evaluated  along 
the  deformed  contour.  Hence  further  numerical  work  is  necessary 
to  quantify  the  effect  of  convective  instabilities  on  far-field 
sound  pressure  and  acoustic  intensity.  In  addition,  care  is 
required  in  attaching  any  physical  significance  to  the  intensity 
vector  plots  (Figures  12  to  15)  at  combinations  of  flow  speed  and 
frequency  giving  rise  to  a  convective  instability. 

According  to  Brazier-Smith  &  Scott  [19],  the  causal  approach 
to  stability  predicts  a  critical  flow  speed  above  which  the  line- 
excited  plate  wiLh  flow  is  absolutely  unstable.  The  flow  speed  is 
given  by 


U  =  0-074pV(D/Ps3ha) 

which  for  a  steel  plate  in  water  reduces  to  a  flow  speed  of 
approximately  15  m/s  (30  knots*)  or  M=0’01.  This  speed  can  be 
determined  by  obtaining  the  dispersion  relation  Z(cj,cr)=0  for 
incompressible  flow  and  considering  when 


az/a<*  =  a2z/aa2  =  o 


i.e.  it  is  equivalent  to  the  case  where  there  are  three 
coalescing  poles  on  the  real  a-axis.  One  consequence  of  this 
onset  of  absolute  instability  is  that  at  flow  speeds  above  the 
critical  speed,  the  pressure  and  acoustic  particle  velocities 
cannot  be  determined  using  the  current  linear  theory.  It  is  well 
known  that  in  a  finite  system  wiwh  flow,  if  a  resonant  mode  of 
the  system  is  excited  at  a  particular  flow  speed  the  level  of 
radiated  sound  increases  dramatically .  This  may  also  be  the  case 
in  situations  where  linear  theory  predicts  absolute  instability 
of  an  infinite  system.  Hence  care  is  required  in  interpreting  any 
numerical  results  at  the  higher  flow  speed. 


9.  CONCLUSIONS  AND  SUGGESTIONS  FOR  FURTHER  RESEARCH 


Expressions  have  been  obtained  for  the  plate  response, 
acoustic  pressure  and  acoustic  intensity  vectors  of  point-  and 
line-excited  infinite  fluid-loaded  plates  in  the  presence  of 
uniform  inviscid  flow.  The  expressions  were  obtained  from  the 
steady-state  equations  of  motion  in  the  form  of  Fourier  integrals 
along  the  real  wavenumber  axis.  Two  expressions  have  been 
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extracted  from  the  literature  for  the  acoustic  intensity  vector 
in  the  presence  of  a  uniform  flow.  For  the  particular  case  of 
steel  plates  in  water  at  flow  speeds  of  less  than  15  m/s  (30 
knots),  the  flow  generally  has  very  little  effect  on  radiated 
sound  pressure  and  the  response  of  a  line-excited  plate  consists 
of  waves  propagating  away  from  the  excitation.  At  very  low 
frequencies,  however,  the  downstream  response  of  a  line-excited 
damped  or  undamped  plate  is  convective ly  unstable  and  one 
component  of  the  upstream  response  is  an  incoming  wave  from 
infinity  which  is  apparently  destabilized  by  damping.  These 
results  are  consistent  with  those  obtained  by  examination  of  the 
dispersion  relation  to  find  free  waves  with  real  wavenumbers 
corresponding  to  complex  values  of  frequency  with  positive 
imaginary  parts.  At  flow  speeds  above  15  m/s,  a  causal  approach 
to  stability  indicates  that  the  coupled  system  is  absolutely 
unstable  and  linear  theory  breaks  down. 

There  are  several  areas  in  which  the  work  can  be  extended. 
Although  the  interesting  physical  effects  occur  at  extremely  low 
frequencies  for  steel  plates  in  water,  for  other  materials,  e.g. 
soft  and  hard  rubber,  the  various  critical  frequencies  and  the 
critical  flow  speed  are  generally  different,  and  the  critical 
frequencies  may  well  be  much  higher,  and  hence  of  more  practical 
significance.  The  theory  can  also  be  extended  to  consider  finite 
plates  or  infinite  strips,  the  stability  of  which  has  been 
tackled  using  the  Benjamin-Landahl  approach  but  not  the  causal 
approach.  There  has  been  much  intramural  and  extramural  work  done 
by  AMTE  and  contractors  on  sound  radiation  and  wave  propagation 
in  more  complex  fluid-loaded  structures  such  as  stiffened  plates, 
infinite  cylinders  with  internal  and  external  fluid  loading,  and 
layered  media.  The  effect  of  flow  on  these  systems  needs  to  be 
studied  and  quantified. 

One  limitation  in  the  present  work  is  that  the  mean  flow  is 
assumed  to  be  uniform.  In  practice  there  is  a  thin  boundary  layer 
next  to  the  plate  surface  in  which  the  mean  flow  velocity  changes 
from  zero  at  the  undisturbed  plate  surface  to  the  free  stream 
value.  In  general  this  boundary  layer  is  turbulent  which  means 
that  the  flow  itself  is  an  additional  excitation  source.  The 
effects  of  laminar  and  turbulent  boundary  layers  on  plates  and 
more  complex  dynamical  systems  need  to  be  examined  and 
quantified.  A  number  of  investigations  have  already  been 
undertaken  on  this  topic.  Benjamin  [14]  has  incorporated  a 
representation  of  a  thin  laminar  boundary  layer  with  zero 
pressure  gradient  in  the  form  of  a  small  correction  to  the 
surface  pressure.  Ff owcs-Williams  &  Purshouse  [22]  have  used  a 
vortex  sheet  model  of  a  boundary  layer  in  the  problem  of  a 
baffled  piston  vibrating  beneath  a  boundary  layer.  They  showed 
that  whereas  a  uniform  flow  at  low  speed  reduces  the  resonant 
frequency  of  the  piston  and  above  a  certain  speed  induces 
instability  and  flutter,  the  vortex  sheet  model  predicts  that  in 
many  cases  the  effect  of  the  flow  is  to  increase  the  resonant 
frequency  without  any  flutter.  These  latter  predictions  agree 
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with  experiment. 


Although  there  has  been  much  theoretical  work  done  on  the 
stability  of  flow  over  surfaces,  there  has  been  relatively  little 
corresponding  laboratory  experimental  work.  The  experiments  of 
Gupta  et  al.  [17]  are  described  in  the  Introduction.  Hansen  et 
al.  [23]  have  conducted  stability  experiments  on  a  finite  plate 
and  rotating  disc.  They  observed  three  different  types  of  surface 
waves  on  the  plate  and  two  on  the  disc.  One  of  the  plate  waves 
was  attributed  to  an  irregularity  in  the  plate  geometry,  but  all 
the  other  waves  could  be  predicted  theoretically.  The  propagation 
velocities  for  all  the  disturbances  were  substantially  below  the 
free-stream  fluid  velocity. 

Because  of  the  low  wavenumbers  and  hence  large  wavelengths 
involved  even  with  very  thin  steel  plates,  experiments  to  test 
the  present  theory  using  steel  plates  in  water  would  be 
impractical.  However  there  may  be  other  materials  with 
substantially  different  physical  properties  which  could  be  used 
in  water -tunnel  tests. 


D. J. Atkins  (HSO) 

Manuscript  completed  August  1982 
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The  general  expression  for  the  acoustic  pressure  amplitude 
P(x,y,z )  of  a  point-excited  plate  is  given  as  equation  (8) 


P(x,y, z)  -  (-ipw^Fo/47r2) 


Z  (vQ,a,0) 


dadj3  (Al) 


where  Z  is  defined  by  equation  (10),  as 


Z(u  ,a,0)  -  D(a2  +  02)2  -  pahu2  -  ipw2 ( 1-aM/k) 2 /y 

0  O  0  0 


with 


y  =  Vt  (k-oM)  2-a2-02) 


In  this  Appendix,  a  non-standard  spherical  co-ordinate  system 
defined  by  equation  (32)  and  illustrated  as  Figure  2(a)  is  used, 
viz : 


x  =  R'cose,  y  »  R’Sin©’Sin^,  z  *  R*sine*cosir 


The  double  integral  given  by  equation  (Al)  earn  be  reduced  to  the 
form 


0)  exp[  iR-g(a,/3)  ]  dad/3 


(A5) 


According  to  the  method  of  stationary  phase,  as  R-*»,  then 


P  «  27ris -f  (a0  ,£0) -exp[  iR-g(a0 ,0O)  ]/[RVI  ab-d2  I  ] 


(A6) 


where  the  stationary  phase  point  (a0,/80)  is  given  by  the 
simultaneous  solution  of 


3g/3a  -  3g/3£  -  0 
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Substituting  equation  (A15)  into  equation  (A12)  gives 

dg/da  =  cos©  -  [ {a+(k-aM)M}/k' ] -sin©  (A16 

which  results  in 

tame  =  V[ (k-a0M) 2-a0  2 ]/[ a0+ (k-a0M)M]  (A17 

This  equation  has  a  closed-form  solution 

a0  =  [k/(l-M2)][-M  +  cos  ©  /  V( 1-M2s in2©) ]  (A18 

Substituting  equation  (A18)  into  equations  (All)  and  (A15)  gives 


©0  -  k •  s in© •  s in?  /  V(l-M2sin2©)  (A19 

In  addition, 

v(aQr0o)  *  k-sine-cosir  /  V(l-M2sin2©)  (A20 

and 

g(aa  rj30)  =  [k/(l-M2) J [-M-cos©  +  V(  1-M2s in2©)  ]  (A21 

The  second  derivatives  of  g  are  given  by 

aag/0a2  -  -sin© • cos?  [k2-( 1-M2 )©2 ]  /  (k'2-©2)3/2  (A22 

a2g/3©2  -  -sin©-cosy-k' 2/(k' 2-©2)3/2  (A23 

a2g/aa3©  -  -sin©-cosy*©[a+(k-aM)M3  /  (k'2-©2)3/2  (A24 


By  evaluating  the  expressions  (A 22)  to  (A24)  at  (aa,£0), 
it  can  be  shown  that 


ab-d2  =  (1-M2sin20)  /  (k2s in20 • cos2**) 


(A25) 


b  -  -  V(l~M2sin20)  /  (k*cos2if) 


(A26) 


Prom  equations  (A26)  and  (A27),  it  is  evident  that  ab-d2>0  and 
b<0.  Hence  a<0  and  s=*-l,  from  which  equation  (A6)  becomes 


P  «  (-2rrik)  •  sin© •  cosy*  exp[  iR* g(k,M,  0)  ]•  f  (a0  ,/30  )/( 1-M2sin20)  (A27) 


By  expressing  f(a0,£0)  and  g(a0»£0)  in  terms  of  k,  M,  0  anc 
f,  some  algebraic  manipulation  gives  the  result  (equation  (36)) 
for  the  far-field  pressure. 


Traditional  Spherical  Co-ordinates 

The  co-ordinate  system  is  defined  by 
x  -  R- sine -cos*,  y  -  R*  s ine • s in*,  2  -  R'cosfl 
The  stationary  phase  point  (ao,j0o)  is  given  by 

a0  -  [k/( 1-M2 ) ] [-M  +  sine-cos*  /  V{h(M,e,*)}] 

-  k-sine-sin*  /  V{h(M,e,*)} 

where 

h(M,e,*)  -  1  -  M2  ( 1  -  sin2e-cos2*) 

The  far-field  pressure  is  given  by 


P(R, e ,*) 


Fofj^l<I-%M/lt)Mh7§T97?i'l9xp{iR,g(1<:-M-8-1,'>1 


(l-aflM/k) 


2  r  ik-cose  1  r  .  ,  2  .  _2  *  2  ..A  , 
L- V{h(M~,e,*)}J {1  <*o*V  /kp} 


where 


g(k,M,e,*)  -  [k/(l-M2) ] [-M- sine • cos*  +  V{h(M,e,*)}] 
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